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ABSTRACT 

We make use of publicly available results from N-body Millennium Simulation to create 
mock samples of lensed supernovae type la and core-collapse. Simulating galaxy-galaxy 
lensing we derive the rates of lensed supernovae and find than at redshifts higher that 
0.5 about 0.06 per cent of supernovae will be lensed by a factor two or more. Future 
wide field surveys like Gaia or LSST should be able to detect lensed supernovae in 
their unbiased sky monitoring. Gaia (from 2013) will detect at least 2 cases whereas 
LSST (from 2018) will see more than 500 a year. Large number of future lensed 
supernovae will allow to verify results of cosmological simulations. The strong galaxy- 
galaxy lensing gives an opportunity to reach high-redshift supernovae type la and 
extend the Hubble diagram sample. 
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1 INTRODUCTION 

Gravitational lensing is a successful tool allowing observing 
surveys to reach to potentially smaller, fainter, and less mas- 
sive objects as it does not rely on extremely deep imaging. 
In strong lensing a massive foreground galaxy deflects the 
light from a background object, resulting in multiple im- 
ages of the source being seen. In additi on the source is typ i- 
cally magnified by a factor of ~10 (e.g. lNewton et ai1l201lft . 
This fact allows for studying of obj ects which are fainter and 
small e r than otherwise possible |Schneider. Ehlers fc Falcq 
1 19921 : ISchneiderl 12006ft . We currently know a few hun- 
dred of strong lens systems. Many such s ystem have 
been discovered in the radio-domain ( e.g. IWinn et all 
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Supernovae (SNe) are divided into two main types I and 
II based on their spectra, with a secondary distinction based 
on light-curve shape. Supernovae type la (SNe la) are widely 
used as a cosmological tool to analyse th e Hubble diagram 
and to constrain parameters Om,^a , w |Riess et al 1 1 19981 : 
lAstier et all 120061 : IConlev et all l201ll : ISullivan et aljkoilft . 



l200d : iBrowne et all 120031 : iMvers et all 12003ft . through 
spec troscopic or colour selection from differe nt surveys. 



e.g. IBolton et all ll2006ft. iBolton et all (|2008r i. iTreu et al I 



d2006ft (SLAC S). iBelokur ov et all (I2009II (CASSOW ARY) 



Brownstein et alj J2012ft (B ELLS). iLimousin et all (|2009ft , 
Gavazzi et all (12012ft (SL2S), lTreu et all |201lft (SWELLS). 



Negrello et al] l|2010ft (Herschel AT L A S ) and in Hu bble 
Space Telescope COSMOS survey (|Faure et all l200Sft at 
optical wavelengths. We observe quasars, radio sources 
and various types of galaxies as source population, which 
provide a wide range of applications of the lensing phe- 
nomenon. Strong lenses have been used, for example, to 
investigate the density profile of lensing galaxies, the evo- 
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SNe la remain, at present, one of the most direct and mature 
method of probing the dark energy. Thought to be the result 
of the thermonuclear explosion of a white dwarf star exceed- 
ing the Chandrasekhar mass limit by accreting mass from 
its companion star in a binary system, or by merging with 
another white dwarf, they are standardisable candles which 
explode with almost the same brightness everywhere in the 
Universe due to the uniformity of the triggering mass and 
hence the available nuclear fuel. Their cosmological use takes 
advantage of simple empirical relations between their lumi- 
nosity and other parameters. The application of relations 
between SN la light-curve-shape, colour, and host galaxy 
properties provides robust distance estimates which allow 
SNe la to be used to measure cosmological parameters. The 
SNe la samples are now sufficien tly numerous that w e are 
able to measure dark energy (e.g. ISullivan et al.ll201lft 

The core-collapse supernovae (SNe CC) are thought to 
be the final evolutionary phase of massive stars with ini- 
tial mass larger than 8 Mq. The progenitors of CC super- 
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novae are massive stars, either single or in binary systems, 
that complete exothermic nuclear burning, up to the devel- 
opment of an iron core that cannot be supported by any 
further nuclear fusion reactions or by electron degeneracy 
pressure. Due to the brief lifetime of their progenitor stars 
(shorter than the typical cosmological time scale), the rate 
of SNe CC occurrence closely follo ws the current star for - 
mation rate in a stellar system (e.g. iBotticella et al.l l2012). 
SNe CC are responsible for a significant fraction of heavy el- 
ements present in the Universe and their energetic impacts 
perhaps trigger further star formation. They are also rele- 
vant to many astrophysical issues, for example, associated 
with the physical and chemical evolution of the galaxies and 
the local environment, the production of neutrinos, cosmic 
rays and gravitational waves. Moreover, SNe CC seem to 
be particularly promising to measure cos mological distance, 
in add ition to type la SNe according to iPumo fc Zampieril 
<|201ll ). 

Forasmuch as the light collection capabilities of tele- 
scopes are limited, our prospect to study high-redshift su- 
pernovae is strongly tentative. The high-redshift SNe la are 
necessary to explore cosmological parameters more accu- 
rately. The strong lensing acting as a cosmic telescope can 
provide an extra insight into higher red shift supernovae and 
open a new area of thi s research (e.g. iGoobar et al.l [20021 ; 
lOguri fc Marshallll2010l ). 

In this paper, we present predictions for expected num- 
bers of supernovae type la and core-collapse which are grav- 
itationally lensed in the most common galaxy-galaxy lens- 
ing. We investigate how many lensed supernovae might be 
detected in currently on-going and future large-scale surveys 
like Gaia or LSST. 

For many of the existing and planned surveys explor- 
ing the transient sky detecting supernovae is one of the main 
scientific goals. The Supernovae Search in Stripe 82 of SDSS 
programme yielded about 500 detections of Supernovae and 
contains almost 4 million stellar objec ts and galaxies mon- 
itored in the it, g, r, i, and z bands (Bra mich et al.l [20081 ; 
lHan et aljEffllcl ). 

ESA's Gaia mission will collect astrometric, photomet- 
ric, spectro-photometric and spectroscopic measurements 
for one billion celestial objects from the entire sky down to 
V ~ 20 mag. The mission is designed to last 5 years, with a 
possible one- year ext ension, after a launch currently p lanned 
for September 2013 (| Wyrzykowski fc Hodgkir]|2012i ). Gaia 
is expected to detect about 6300 SNe (to 19 mag), of which 
85 per cent will be las, of which 30 per cent will be detected 
before maximum jBelokurov fc Evans! [20031 ; I Altavilla et al.l 
120121 ). The maximum redshift detectable by Gaia for SN 
type la is about 0.15, however availability of low-resolution 
spectra for all objects will significantly decrease the rate of 
false-alerts and will all ow for immediate superno vae typing 
and redshift estimation (|Wvrzvkowski et al.l20 12.; Blagorod- 
nova et al. in prep.). 

The most ambitious large-scale survey currently 
planned in the optical domain is the La rge Synoptic Sur- 
vey Telescope (LSST) (|Abell et al.ll2009l ). LSST will be a 
wide-field ground-based system designed to obtain multiple 
filter images covering the sky visible from Northern Chile. 
The current baseline design, with an 8.4m (6.7m effective) 
primary mirror, a 9.6 deg 2 field of view and a 3.2 Gigapixel 
camera, will allow for about 20,000 square degrees of sky in 



total to be covered. One of the ma in scientific goals of LSST 
is exploring transient optical sky (jlvezic et al.ll2008l ). 



This paper is organised as follows. In Section 2 we de- 
scribe the Millennium Simulation with its applications. We 
then present supernovae rates used in our simulations in 
Section 3 and in Section 4 we describe lensed supernovae 
simulation and show our results. After that we discuss the 
results and conclude in Section 5. Throughout this paper 
we assume a flat A-cold-dark-matter (ACDM) cosmologi- 
cal model with parameters J1a = 0.73, Om = 0.27 and 
Ho = 71km s _1 Mpc _1 , h = 0.71 as favoured by the seven 
year Wilkinson Microwave Anisotropy Probe (WMAP) re- 
sults |Komatsu et al.|[201ll 'l. 



2 MILLENNIUM SIMULATION AND ITS 
APPLICATIONS 

In this study we rely on the Millennium Simulation 
ijSpringel et al.|[2005h Fl which offers high spatial and time 
resolution within a large cosmological volume, ideal for 
studying lensing at various scales. The Simulation describes 
the evolution of matter distribution and formation of the 
cosmic structures in a cube of co- moving size 500ft~ 1 Mpc 
and extends the possibility of realistic modelling of the 
light propagation in the Universe. In particular the mat- 
ter distribution is described by the positions of the ~ 10 10 
simulation particles known at 64 epochs representing red- 
shift from to 127. The positions, virial masses, virial 
radii, and virial velocities of gravitationally bound haloes 
(containing at least 20 s i mulation partic l es) o btained by 
|Pe Lucia fc BlaizotJ (120071 ). iBertone et all |2007l) are taken 
from the Virgo - Millennium Database website^ In the whole 
study we consi der the haloes as the s i ngular isothermal el- 
lipsoi ds (SIE) (|Kormann et al.l 1 1994 iKeeton fc Kochanekl 
Il997h . Every halo from the simulation is a probable lens 
which will be described by a few parameters: the Einstein 
radius, the axis ratio, the position angle. The Einstein ra- 
dius depends on the mass of the dark matter halo (or here 
on the virial velocity of the halo), the distance between the 
observer and the source and the distance between the lens 
and the source: 



Osi 



4ty 



cr 2 D ls 



(1) 



1 In the Millennium Simulation the ACDM cosmological model 
with parameters obtained by WMAP 1-year-observations (Qa = 
0.75, f?M = 0.25 and Ho = 73 km s~ 1 Mpc~ 1 ) was assumed. For 
this kind of simulation (specifically the analysis of gravitational 
instability) two parameters play an important role: the amplitude 
of density fluctuations, erg = 0.9 (as compared with the current 
value of 0.8), and the density perturbation spectral index n = 

1 (0.963 respectively). The higher value of parameter erg in the 
simulation implies that the gravitationally bound structures in 
the mechanical equilibrium appear earlier there than they did in 
reality. The discrepancies in the evolution between the simulation 
and the Universe increase with the object mass. We are not able to 
transform the results from the simulation to comply with todays 
knowledge. We are aware that the lensing probability may be 
overestimated, although for supernovae bursting at redshifts of 
less than 2.5 (like in our study) the effect should not be dramatic. 

2 http:/ /gavo. mpa-garching.mpg.de/Millennium/ 
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Figure 1. Upper panel: The normalized distribution of the virial 
velocity of dark matter haloes in the Millennium Simulation. 
Lower panel: The contribution of each velocity bin to lensing 
probability. Shaded region - the range of typical galaxy virial 
velocities. 

where a is the velocity dispersion of the halo, c is the light 
speed, q is the axis ratio, Di 3 ,D oa are cosmological angular 
diameter distances from the lens to the source and from 
the observer to the source, respectively. The haloes in the 
simulation are described by the virial velocity which satisfies 

the relation for rotation velocity: v V i r = y^p^S where 
M V i r is the virial mass, R v i r is the virial radius, G is the 
gravitational constant. Hence, by taking into account the 
relation v 2 otation — 2a 2 , the Einstein radius in our study 
can be expressed by: 

e S iE = 2^%^. (2) 

We, therefore, use the results from the Millennium Sim- 
ulation database to obtain the distributions of feasible lens 
properties and the lensing probability. 

2.1 The lens velocity distribution 

For each redshift epoch we obtain the virial velocity distri- 
bution of gravitationally bound haloes versus redshift which 
we use to randomise the Einstein radius. As it can be seen 
in Figure [T] (upper panel) that the velocities of the haloes 
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Figure 2. The whole-sky averaged galaxy-galaxy lensing proba- 
bility as a function of source redshift. 

extend from 20 to 1000 km s _1 . In this study we only investi- 
gate lensing by individual galaxies, as it is far more common 
than galaxy cluster lensing. Indeed, there are many advan- 
tages of observing massive galaxy clusters to find lensed su- 
pernovae due to their large regions of high magnification. 
These asp ects were already fu lly discussed in lGoobar et al.l 
(2009) and lRiehm et al.l (|201ll ). We assume that the velocity 
dispersion of a typical galaxy should be in range between 70 
and 400 km s" 1 (100 < v vir < 570 km s" 1 ), following the 
fitting o f a modified Schec hter function to the observed SDSS 
sample l|Choi et al-l feOQ?), This restricts us from using very 
massive haloes (which can be considered as galaxy clusters) . 
Moreover, we neglect a group of low-mass haloes (expected 
as dwarf galaxies), as on these scales (iw < 100 km s _1 ) 
cold dark matter simulations appear to overestimate the 
amount of structure, to a de gree that appears in consistent 
with current observations (e.g. lSpringel et alj2008i ). We can- 
not accurately take into account the effect of dwarf haloes on 
lensing probability, because there is insufficient data on their 
real distribution. We therefore exclude them altogether from 
this study. However, one can estimate that their inclusion 
would increase the lensing probability by 6-10 per cent. 

2.2 Galaxy-galaxy lensing probability 

In con structing the matter layers we follow ICarbone et al.l 
(2008), randomly rotating and shifting simulation cubes cor- 
responding to different epochs, which eliminates the conse- 
quences of periodic boundary conditions used in their cal- 
culations. If we consider a source at epoch n, then there are 
(ri — 1) layers between the observer and the source, where 
the lens can be situated. We select ten independent direc- 
tions on the sky and examine fields covering 3° x 3° each. 
Each halo from the cone contributes to the strong lensing 
probability as 7rO| IE /(3° x 3°), see the lower panel of Figure 
fT] Afterwards, we average the probabilities from each of the 
ten fields to estimate the probability of strong lensing on the 
whole sky. As a result we obtain the strong lensing probabil- 
ity for each source epoch, represented by the epoch redshift, 
see Figure [5] Lensing probability rises sharply in the red- 
shift regime from to 0.5, reaching 10 -5 and then flattens 
at about 10 -4 to 10 -3 for higher redshifts. We acknowledge 
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Figure 3. The lens location probability in the Millennium Sim- 
ulation for sources at different rcdshifts. 



Figure 4. The shear distribution in the Millennium Simulation 
for the sources at different redshift. 



the fact that by treating each halo as a SIE we add extra 
mass to the line of sight. That leads to an overestimation of 
the Universe density and therefore of the magnification and 
lensing rates. We examined a few hundred haloes randomly 
chosen from the Millennium Simulation compensated by a 
constant, negative surface matter density disk (correspond- 
ing to the convergence k) and we claim that the cross-section 
for multiple images does not change in such a model. The 
separation of images changes by a factor 1/(1 + |k|) and the 
amplification by 1/(1 + | ac | ) 2 , where the value of parameter 
k| is really small (< 10 -3 ), hence we conclude the problem 
of overdensity can be neglected in this study. 



2.3 The lens location probability 

After we choose the source redshift, we are interested in the 
likelihood of a lens presence in every redshift bin between 
the observer and the source. For every source epoch we saved 
the contribution of every redshift layer to the total proba- 
bility, so we are able to randomise the lens redshift with this 
distribution. Figure [3] shows the lens location likelihood for 
sources at three different redshifts. 



2.4 The shear distribution 

The smoothed distribution of matter and the positions and 
properties of gravitationally bound haloes are sufficient to 
describe the light propag a tion i n terms of convergence. In 
IJaroszvriski fc Kostrzewal (|2010l ) use of the Born approxi- 
mation enabled the construction of the mass distribution 
and the calculation of the convergence and shear maps. We 
utilise those results to obtain the distribution of the shear, 
which will contribute to the lensing signal as an effect of the 
lens environment (Figure 



3 SUPERNOVAE RATES 

The volumetric supernova type la rate at the cosmic time 
t - Ria(t) is the convolution of the cosmic star formation 
history (SFH) S(t) with the delay-time distribution (DTD) 
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Figure 5. The rates of non-lensed super novae type la for different 
SFH function ffrom lGraur et all ll201ll) '). 
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Figure 6. The ra tes of non-lensed C C supernovae for SFH func- 
tion Y only ffrom lGraur et al.l <201ll) '). 
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Table 1. Star formation histories and results of best-fitting DTD 
parameters. 
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Figure 7. The rates of lensed (lower curves) and non-lensed (up- 
per curves) supernovae of SNe la (solid lines) and SNe CC (dashed 
lines). The rates were calculated supposing the SFH Y function. 



*(*)= 
Ru(t) 



S(t - r)*(r)dr. 



(3) 



The SFH function has been measured out to re dshift 8 
i jYuksel et al.ll200ct l. Following [Graur et al l l|201ll ) we use 
four different SFHs in order to investigate any dependance 
of lensed supernovae rates on the SFH function. The SFH 
functions are defined as broken power laws in the form of 
S(z) ~ (1 + z) 1 , as propo sed bvlYiiksel et alj (|2008l ) (here- 
after denoted as Y) and lOda et al.l ( 20081 ) (hereafter Ol 
and Q2) and S(z) ~ (a + bz)/(l + (z/c) d ), as defin ed by 
iHopkins fc Beacoml |2006l ) (HB). iGraur et alj l|201ll ) show 
that only power-laws of the form ^(t) ~ (t/lGyr) -13 for the 
DTD function are consistent with the da ta. Table [J contain s 
best-fitted values of all parameters from I Graur et all (|201ll ) 
for all SFH definitions used here. Different SN la rates re- 
sulting from different SFHs are presented in Figure [5] 

For the core-collapse supernovae we also use a rescaled 
SFH function because of the suggestion that this rate is sim- 
ply proportional to the cosmic star formation rate. In Figure 
[6] we plot the core - collap se SN rate with the normalisation 
from IGraur et alj |201ll ). Moreover, in Figure [7| we show 
the comparison between the rates of lensed and non-lensed 
supernovae of both types. 



4 LENSED SUPERNOVAE 

We choose the singular isothermal ellipsoid (SIE) model with 
external shear as the lens model for haloes from the Millen- 
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Figure 8. Upper panel: The magnification distribution for 
sources at redshifts 0.5, 1.5, 2.5. Lower panel: The magnifica- 
tion versus lens virial velocity (used for estimating the mass) for 
sources located at redshifts 0.5, 1.5, 2.5. 
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Figure 9. The cumulative number of lensed SNe la per year per 
whole sky (peak magnitudes in rest-frame B-band) for different 
SFH. 
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Figure 10. The cumulative number of lensed SNe CC per year 
per whole sky (peak magnitudes in rest-frame B-band). 
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nium Simulation. The SIE model is parametrized by the 
angular Einstein radius Osie, which is related to the virial 
velocity v V ir of the lens galaxy that we randomise based on 
the distribution obtained from the Millennium Simulation. 
We assume a Gaussian distribution for th e axis ratio with a 
mean of 0.7 and dispersion of 0.16 as in lOguri fe Marshalll 
The distribution is truncated at q = 0.5 and q = 1. 
Using simulated data we model about 350,000 lensing 
event s of SNe sources in the GRAVLENS programme IjKeetonl 
l200ll ). For each system we obtain the image positions and 
the image magnifications of sources randomly situated in the 
background of lenses. The position of the source (x,y rela- 
tive to the lens axis) fulfils the condition x 2 + y 2 < Os IE . By 
choosing the image with the highest magnification we can 
calculate the brightest observed magnitude of this lensed su- 
pernova. It should be stressed that no a priori assumptions 
were made about the resulting images e.g. magnification, 
separation etc. However, because the source is located inside 
the Einstein radius, there is more than one lensed image in 
the strong lensing case. Our research concerned the lensing 
statistic, hence we are content with detecting even just one 
image at a given magnitude threshold. From this simulation 
we obtain the amplification distribution for sources at dif- 
ferent redshifts, see upper panel in Figure [8] As expected, 
the most probable magnifications are near 2, but there is 
a non-negligible probability of obtaining magnifications in 
excess of 10 or 100. The highest magnifications are due to 
specific source-lens system parameters, especially the posi- 
tion of the source, but not due to the source redshift, what 
was investigated in three separate redshift bins, shown in 
Figure [H] Lower panel of this figure shows that a source can 
be highly magnified even by a dark matter halo with non- 
high velocity (and therefore mass), which can occur at all 
redshifts. 

Supernovae type la are thought to have the same value 
of absolute magnitude Mb at the maximum of the light 
curve in B-band. This fact is ascribed to the origin of these 
objects as we mentioned in introduction. However, the same 
cannot be said of core-collapse supernovae. Hence, we as- 
sume the average absolut e magnitude of a supernova type 
la as M B = -19.13 mag l|Conlev et al.ll201ll ; ISullivan et all 
[loll]). 

Core-collapse SNe have a wide range of absolute mag- 
nitudes and we choose Mb = —18.5 mag, as an average 
value base d on SN light curve t emplates from iNugent et all 
|2002l ) and lBarbarv et al l (|2012l ). The lensed supernovae will 
be magnified by a factor —2.5 log(/i) where in our case /x is 
the largest magnification in the simulated lens system. The 
peak of the observed light curve peak has rest-frame B-band 
magnitude of 

mT k = Mb + 51og _ 2 . 51og(At)) (4) 

where di,{z) is the luminosity distance. 

Combining the SN rates and the magnification distri- 
bution we obtain lensed supernovae rates for both types: 
la and core-collapse. First of all, Figure [5] shows that we 
expect about 0.6 (between 0.53 and 0.66 for different SFH 
functions) lensed supernova type la per year up to magni- 
tude 20 on the whole sky. However, for future surveys (with 
24 magnitude capability) we predict about 700 detectable 
SNe type la per year on the whole sky. 



Nonetheless, a closer look at the rates histogram reveals 
that those measurements cannot be expected to easily dis- 
tinguish between different SFH functions. As seen in Fig. 
SFH functions differ significantly only at higher redshifts, 
i.e. above 1.5. As most SNe are magnified only by a small 
factor these redshifts correspond mostly to SNe of magni- 
tude above 24.5, which is near the limit of our study. 

For the core-collapse SNe we predict 1.3 and less than 
1700 lensed objects per year for the whole sky within 20 and 
24 mag range, respectively (Figure [TP)). 

In this study we focus on detectability of the strongest 
magnified image of a supernova. However, when the separa- 
tion between the lensed images is smaller than the angular 
resolution of a detector, the magnified images will merge, 
forming a single blended observed image. We tried to es- 
timate the brightness of such blended images taking into 
account both temporal evolution of the light curve of the 
lensed supernova and the time delay between the images. 
In our sample of lenses most image separations fall below 
0.5 arc seconds for double-imaging lenses and below 0.2 arc 
seconds for quads, which is comparable with the resolving 
powers of typical ground-based surveys. For the two-image 
systems (more than 97 per cent of systems in whole sample), 
only in case of 8 per cent the shift in observable brightness is 
significantly higher than in case of a single image by at least 
of 0.5 mag. The average shift in brightness is 0.23±0.17 mag, 
which is comparable with the error in the estimation of the 
absolute magnitude for the SNe of both types. In case of the 
four-image systems the overall brightness of blended images 
can be much larger than the brightness of the strongest im- 
age. The average shift in brightness is about 1 mag, reaching 
up to 1.4 mag, however, the quads compose only 3 per cent 
of our sample, hence the overall effect of blending of the 
images can be neglected. 

We do not take into account microlen sing of lensed SNe 
cause d by the stars of the lens galaxy (|Dobler fc Keetorj 
2006) as this would be beyond the scope of this work. We 
also neglect the effect of extinction in the foreground galaxy, 
which is claimed to affect neither t he detectability of l ensed 
events nor the flux measurements dGoobar et alj|2002h . 

Owing to the fact that we are interested in the statis- 
tical probability of lensed supernovae appearance we do not 
investigate the impact of the source population characteris- 
tic. Obviously this seems to be important during choosing 
specific lens systems for observation. 



5 DISCUSSION & CONCLUSIONS 

In this paper we present the simulations of supernovae be- 
ing lensed by galaxies, represented by haloes from Millen- 
nium Simulation. The numbers of expected lensed super- 
novae in large scale surveys are gathered in Table[2]for SDSS 
Stripe82, Gaia and LSST. Given their limiting magnitude, 
assuming about 20 per cent detection efficiency we estimate 
that in the entire Stripe 82 data (about 8 years) there is 
hardly any chance (10 and 20 per cent for SN la and CC, 
respectively) that the data contains any lensed supernova, 
unless we got really lucky. 

Gaia mission will cover the entire sky in its 5 years 
(2013-2018) regular observations. Its sampling will depend 
on the region of the sky, hence the efficiency will vary. 
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Table 2. The expected number of lensed SNe type la and core-collapse SNe in various surveys (expected efficiency about fOO and less 
than 20 per cent). 



Survey 


Survey area 


Survey time 


Survey limiting 


Number SNIa 


Number SNCC 


Number SNIa 


Number SNCC 




[dcg 2 ] 


[years] 


magnitude (filter) 


[1 year] 


[1 year] 


[survey time] 


[survey time] 


Stripe 82 (100%) 


300 


8 


21.5 (i) 


0.073 


0.124 


0.584 


0.992 


Gaia (100%) 


40000 


5 


20.0 (V) 


0.60 


1.27 


3.0 


6.35 


LSST (100%) 


20000 


10 


24.5 (i) 


710. 


1900. 


7100. 


19000. 


Stripe 82 (20%) 


300 


8 


21.5 (i) 


0.015 


0.025 


0.119 


0.198 


Gaia (20%) 


40000 


5 


20.0 (V) 


0.12 


0.254 


0.6 


1.27 


LSST (20%) 


20000 


10 


24.5 (i) 


142. 


380. 


1420. 


3800. 




-6 



0.0 0.5 1.0 1.5 2.0 

Redshift 

Figure 11. The Hubble diagram for lensed SNe. The solid line represents dist ance modulus for c osmological parameters: Qm = 0.27, 
SIa = 0.73, dotted line denotes the scatter in current combined sample of SNe la llConlev et al.|[201 ll) . boxes - simulated lensed SNe la for 
1-year-observations on the whole sky with 100 per cent efficiency. Lower panel - residuals from the cosmological model. Inset: distribution 
of number of lensed and non-lensed SNe la in the sample at redshift of about 1.5. The solid line repre sents the non-Iensed sample, the 
dashed line - the lensed sample. The arrow points to location of a probable lensed supernova SN1997ff iLewis fc Ibata 2 00 1[ ). 



Assuming 20 per cent average detection efficiency we es- 
timate that during the entire mission there might be 
one lensed supernova detected of each type. Gaia data 
processing pipeline will have the special ability to clas- 
sify supernovae into types and to estimate the redshifts 
based on the low-resolution spec t rograph available on-board 
jWvrzvkowski fc Hodgkml |2012| . IWvrzvkowski et all |2012| . 
Blagorodnova et al. in prep.). Therefore, a lensed case of su- 



pernova should be recognised promptly and will be alerted 
to the community immediately. 

Because of its much deeper imaging capabilities, the 
LSST (planned to operate from 2018-2020) will give us more 
rewarding results with a significant number of detections of 
lensed supernovae. Again, assuming 20 per cent efficiency 
for the part of the sky LSST will cover in 10 years of ex- 
pected operation there should be about 1400 lensed super- 
novae type la and more than 3800 lensed core-collapse su- 
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pernovae. This will significantly enhance our probing of the 
cosmological distance scale with large statistical samples of 
high redshift supernovae. The estimated number of SN de- 
tectio ns are higher than the estimates of lOguri fc Marshall 
(2010) on account of no assumptions being made on detec- 
tion of a second image (e.g. separation and magnitude limit), 
which being essential in that study. 

Our simulation indicates that sources at redshift larger 
than 0.5 have lensing probability better than 10 _J (Figure[2]) 
and vast majority of lensed supernovae originate from these 
redshifts. The rates of lensed supernovae of redshift 0.5 and 
higher are about 10 -9 to 10 -7 , indicating these are still rel- 
atively rare events, even in high redshifts. When compared 
to the rates of all supernovae we expect that at redshifts 
0.5 and higher about 0.06 per cent of all supernovae will be 
lensed at least by a factor 2 (Figure [7J. The majority of the 
brightest images of lensed SNe will be magnified by factor 2 
(i.e. by 0.75 mag), however not negligible part of them will 
have magnifications exceeding 10, making the lensed super- 
novae easily standing out of the bulk of typical supernovae. 

The hunt for lensed supernovae can be also performed 
in a targeted regular observations of selected known galaxy- 
galaxy lenses. Our study demonstrate that the most impor- 
tant factor boosting the possibility of higher magnification 
is the position of the lensed galaxy on the source plane. If 
the supernova host galaxy intersects caustics we will obtain 
higher probability of larger magnification, especially if we 
expect more supernovae to burst closer to the centre of the 
source galaxy. 

Moreover, when looking for supernovae lensed from red- 
shift 0.5, their B-band peak at the rest-frame due to cosmo- 
logical shifting will occur in near-infra red. Therefore, ob- 
serving at these wavelengths will assure the supernova is 
observed at its brightest, hence improving chances of its de- 
tection. 

We also investigate the impact of lensed supernovae 
type la on the Hubble diagram (Fig. Ill|) . Generally, it can 
be claimed that the majority of lensed supernovae sample 
will be distinguishable from the current sample. However, 
every outlier should be checked against being contaminated 
by a lensing event. The inset in Figure fTTI shows whether two 
samples - non-lensed and lensed supernovae type la which 
burst in redshift bin between 1.4 and 1.5 - are distinguish- 
able. Here we as sume average SN m agnitude as 26.1 with 
dispersion ~0.2 (iConlev et al.ll201lh . The dispersion is the 
result of an intrinsic scatter in observations due to an incom- 
plete understanding of SN physics or all photometric uncer- 
tainties. It can be expected that a part of a lensed sample is 
mixed with an non-lensed sample. For example, SN 1997ff 
discovered in the Hubble Deep Field was identified as a high- 
redshift supernova m ildly magnified by gravitational lensing 
l|Lewis fc I bata 2001:). Two galaxies lying close to the line- 
of-sight enhanced the magnitude of this SN by ~0.4 mag (as 
indicated in Fig. Ill[) , a change larger than the differences of 
various cosmological models' predictions. 

Our results are in strong agreement with a recent 
study bv lKarpenka et all (|2012l ). where they demonstrate a 
Bayesian statistical methodology for constraining the prop- 
erties of dark matter haloes of foreground galaxies that in- 
tersect the lines-of-sight towards supernovae type la. The 
method was applied to realistic simulated SNIa data, based 
on the real 3-year Supernova Legacy Survey (SNLS3) data 



release (about 160 SNela at redshifts 0.1 < z < 1 up to 
24.5 mag, on 4 deg 2 survey area). Two dark matter halo 
density profile of foreground galaxies models were chosen: a 
truncated singular isothermal sphere and a Navarro-Frenk- 
White profile. They show that detection of a lensing signal 
in the simulated sample is highly unlikely for both dark mat- 
ter halo profiles. Albeit, in the analysis of real SNLS3 data 
they obtain a very minor detection of the lensing signal (for 
the truncated singular isothermal sphere model only). For 
such a small survey area we expect as few as 0.028 lensed 
event per year (assuming 20 per cent detection efficiency). 

Our study shows that even though some surveys will 
detect large samples of lensed supernovae, the differences 
due to different SFH functions start to emerge only at a 
level of about 26 magnitude, not available yet in wide field 
surveys. However, the robust sample of lensed supernovae 
can give answers to a number of questions generated by the 
cosmological simulations results, e.g. about the universality 
of the dark matter density profiles. 
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